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基于 位 势 理论 的 磁场 预测 法 * 
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摘 要 : 为 了 提高 舰 船 的 隐藏 性 ， 需 要 对 舰 船 进行 消 磁 ， 而 磁场 预测 是 舰 船 消 磁 的 核心 步骤。 对 于 任意 舰 
船 ， 本 文 构造 了 包含 它 的 一 个 封闭 曲面 ， 并 利用 位 势 理论 ， 建 立 该 曲面 上 等 效 磁 荷 面 分 布 与 磁场 
传感器 测量 值 之 间 应 该 满足 的 积分 方程 。 离 散 化 求解 后 ， 利 用 该 等 效 磁 荷 面 分 布 对 舰 船 产生 的 磁 
场 进行 预测 。 此 方法 既 不 必 考 虑 被 研究 舰 船 的 几何 结构 ， 又 可 以 在 任意 方向 上 进行 预测 ， 因 此 适 
用 范围 广 。 数值 模拟 表明 了 该 方法 预测 的 精确 性 和 应 用 价值 。 

关键 词 : 舰 船 ， 消 磁 ， 磁场， 预测 ， 位 势 理 论 ， 积 分 方程 


分 类 号 : AMS(2000) 31B10 中 图 分 类 号 : 0175.5; U665.18 文献 标识 码 : A 
1 引言 


现代 舰 船 大 多 是 由 钢铁 材料 建造 的 ， 它 在 地 球 磁场 的 作用 下 ， 犹 如 一 块 巨大 的 浮动 磁体 ， 
在 其 周围 空间 产生 磁场 。 磁 引信 水 雷 、 磁 性 制导 鱼雷 等 磁性 武器 就 是 利用 舰 船 磁场 发 起 攻击 
的 。 因 此 为 了 减 小 舰 船 磁场 ， 提 高 隐蔽 性 ， 需 要 对 舰 船 进行 消 磁 。 舰 船 磁场 预测 问题 就 是 在 舰 
船 周围 放置 若干 个 磁场 传感器 ， 利 用 所 测 得 的 舰 船 磁场 的 强度 信息 来 预测 其 在 空间 中 任意 点 上 
产生 的 磁场 。 此 问题 是 舰 船 消 磁 的 一 个 核心 步 又， 因为 只 有 在 较 精确 地 预测 了 磁场 之 后 ， 才 能 
通过 调整 消 磁 绕组 的 电流 来 产生 与 预测 值 相反 的 磁场 ， 从 而 将 舰 船 保持 在 较 低 的 磁性 水 平 ， 实 
现 舰 船 消 磁 。 常 用 的 磁场 预测 方法 有 深度 换算 法 由 和 磁体 模拟 法 上 3 等。 其 中 深度 换算 法 只 能 
在 半空 间 中 进行 舰 船 磁场 的 预测 ， 而 磁体 模拟 法 需要 提供 舰 船 等 铁 磁 物 体 的 几何 结构 才能 保证 
预测 的 精度 。 文 献 [3] 和 [4] 分 别 用 遗传 算法 和 几何 方法 来 定位 舰 船 所 对 应 的 磁 偶 极 子 ， 但 是 考 
虑 的 都 是 极其 简单 的 情况 ， 离 开 实 际 应 用 还 有 相当 大 的 一 段 距离 。 

本 文 对 任意 的 舰 船 ， 以 Laplace 方 程 Neumann 外 问题 为 模型 ， 利 用 位 势 理论 回 ， 建 立 了 包 
围 该 舰 船 的 封闭 曲面 上 等 效 磁 荷 面 分 布 与 磁场 传感器 测量 值 之 间 应 该 满足 的 积分 方程 。 通 过 离 
散 化 求解 该 积分 方程 ， 得 到 等 效 磁 荷 面 分 布 ， 该 分 布 在 曲面 外 产生 的 磁场 与 被 研究 舰 船 等 铁 磁 
物体 产生 的 磁场 是 相同 的 ， 于 是 便利 用 该 分 布 进行 舰 船 磁场 预测 。 本 方法 统一 考虑 了 舰 船 等 铁 
磁 物 体 的 感应 磁性 和 固定 磁性 四， 可 以 在 任意 方向 上 进行 预测 ， 运 用 时 不 必 提 供 舰 船 等 铁 磁 物 
体 的 几何 结构 ， 因 此 适用 范围 更 广 。 数 值 模拟 的 结果 表明 了 该 方法 的 精确 性 。 


2 位 势 理论 


首先 利用 位 势 理论 ， 给 出 Laplace 方 程 Neumann 外 问题 的 等 价 积分 方程 描述 。 
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定义 Im 在 光滑 的 闭 曲 面 了 上 给 出 连续 函数 J/， 要 求 出 这 样 一 个 函数 yp(z,y,z)， 它 在 工 的 
外 部 区 域内 调和 ， 在 人 QUT 上 连续 ， 在 无 穷 远 处 满足 


Jim g(z,y,2)=0, ro Vs? yt 2 


而 且 它 在 TT 上任 一 点 的 法 向 导数 2e 存在 ， FWE Zr = f. x n EKR 的 外 法 线 ， 即 
曲面 工 的 内 法 线 。 这 样 的 问题 称 为 Laplace 方 程 Newman 外 问题 。 
定义 2/1 enc 


e(M) - z- «| | as, 0 


称 为 单 层 位 势 ， 其 中 7px 是 曲面 上 的 变 点 书 与 空间 中 (包括 曲面 本 身 ) 任意 点 M 之 间 的 距 
离 ，dSp 是 TFT 上 PP 点 处 的 面积 元 。 

单 层 位 势 p(M) 的 物理 意义 就 是 曲面 了 上 的 面 密 度 为 w(P) 的 磁 荷 层 在 空间 中 M 点 处 产生 
的 标量 磁 位 。 

E8155. 若 w(P) 是 光滑 闭 曲 面 了 上 的 分 片 连续 函数 ， 则 

1) 单 层 位 势 (1) 当 点 M 不 在 上 时 处 处 是 调和 的 ， 且 在 无 穷 远 处 正则 

2) "AM eT 时 ， 单 层 位 势 (1) 中 的 积分 仍 是 收敛 的 ， 且 它 在 了 上 连续 ; 

3) rh (1) 所 确定 的 单 层 位 势 p(1M) 在 任 一 点 Py € T 存在 法 向 导数 I, XH 


Op(P) 1 cos(rpp,, np, ) wu (P5) 
np 4r / lj AR PP, A T 
r 


而 且 (2) 式 右 端的 积分 是 收敛 的 ， 其 中 rpp PHA PRUE, rem 是 rpm KE, np 是 
曲面 了 上 及 点 处 的 内 法 线 的 单位 方向 向 量 。 

由 上 述 定理 ， 光 滑 闭 曲面 上 的 分 片 连续 密度 分 布 的 单 层 位 势 在 区 域 Q' 中 确定 了 一 个 调和 
函数 ， 并 在 Q/' UT 上 连续 ， 于 是 可 以 利用 其 法 向 导数 在 TT 上 的 极限 值 ， 来 选取 适当 的 密度 分 
布 w(P)， 使 对 应 的 单 层 位 势 满足 Neumann 外 问题 给 定 的 边界 条 件 

3|. = f(Po) Per. 
考虑 到 等 式 (2)， 得 到 w(P) 所 应 满足 的 积分 方程 为 


(By) — z [f «^ py dSp2/(PR), P ET. (3) 


一 个 弱 奇 性 的 第 二 类 Fredholm 积 分 方程 。 可 以 证 明 积分 方程 (3) 对 于 任意 的 连续 函 
PUE, f), FÆ, Laplace 7; f£ Neumann 外 问题 


Ap(M)=0, Meg, " 
4 
aà£|.-f(F) Py eT 
的 正则 解 可 取 为 
e(M) = > iE zi 2 dS, MeQ, (5) 


EP o(P) 满足 积分 方程 (3)。 
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3 ” 舰 船 磁场 预测 公式 


假设 用 一 个 闭 曲 面 工 把 整个 空间 划分 成 2 和 89' 两 个 区 域 ， 其 中 r 是 工 的 外 部 区 域 。 现 在 要 
预测 处 在 区 域 虽 中 的 舰 船 在 区 域 2 中 产生 的 磁场 ， 见 图 1。 由 于 磁场 等 于 标量 磁 位 的 负 梯 度 ， 
所 以 可 以 把 磁场 的 预测 转化 为 对 标量 磁 位 的 预测 。 男 一 方面 ， 在 实际 工程 应 用 中 ， 目 前 只 有 磁 
场 传感器 ， 而 不 存在 磁 位 传感器 ， 于 是 可 以 把 磁场 传感器 在 上 的 磁场 测量 值 看 作 边界 条 件 ， 
用 Laplace 方 程 Neumann 外 问题 描述 TU W 中 的 标量 磁 位 。 


5b BE DX IR CY 


n s ON 


M cen PT NO 


P / 
/ 
V4 


内 部 区 域 Q 
磁场 传感器 
图 1: 空间 的 划分 
由 电磁 学 知 ， 在 空间 的 每 个 正则 点 上 ( 即 在 该 点 的 邻 域内 介质 的 物理 特性 连续 加 )， 恒 定 磁 


场 基 本 方程 的 微分 形式 为 0V.B = 0，V x H = J， 其 中 BB 为 磁 感 强度 ， 耳 为 磁场 强度 ，J 为 电 
流 密度 。 假 设 区 域 Q 中 没有 电流 分 布 ， 于 是 可 以 引入 标量 磁 位 yp 来 描述 磁场 


H = -Vọ. (6) 
EO 中 的 介质 是 均匀 、 各 向 同性 的 ， 有 
V-B =V- (uH) = uV -H = ~V’ = 0, (7) 
其 中 凡是 介质 的 绝对 磁 导 率 。 因 此 标量 磁 位 vp 在 区 域 Q' 中 满足 Laplace 方 程 
Ay(M)=0, MEW. (8) 
现在 考虑 Neumann 外 问题 边界 条 件 ， 注 意 到 关系 式 (6)， 可 以 得 到 
ael. = Ve. n|, = -H- n'|, = H. n|,, (9) 
其 中 是 曲面 内 法 线 的 单位 方向 向 量 ,，n 是 曲面 外 法 线 的 单位 方向 向 量 。 这 样 就 将 区 
域 Q' 中 的 磁 位 归结 为 一 个 Neumann 外 问题 ， 边 界 条 件 由 闭 曲 面 上 的 磁场 及 TT 的 法 方向 确 
定 


于 是 再 利用 (3)， 在 知道 了 上 的 磁场 强度 旦 (R)，Ro € 后 ， 可 以 构造 上 的 磁 荷 分 
布 w(P)， 它 满足 积分 方程 


1 ; 
w(Po) = —5- // tp) enun) ne) asp +2H(P) np,, Po eT. (10) 
T E TPP 
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这 样 的 磁 荷 分 布 w(P) ERR Y cp P E ft fr ER O 中 的 舰 船 在 区 域 0 中 产生 的 磁 位 是 相 
同 的 ， 因 此 ， 两 者 的 磁场 也 相同 。 得 到 磁 荷 分 布 后 ， 利 用 等 式 (5) 和 (6)， 可 以 求 出 9 中 的 磁 


场 强 度 预测 值 
1 P 
-ven - - s. [[ v(5 
r 


1 TPM ! 
Tx AL M , 11 
二 JI «Ps asp Men (11) 


H(M) 


Il 


3k (10), (11) 即 为 磁场 预测 的 理论 公式 。 

舰 船 磁场 预测 的 理论 公式 (10)，(11) 还 不 能 在 工程 中 直接 使 用 ， 这 是 因为 : 首先 ， 工 程 中 
只 能 在 曲面 荆 的 有 限 点 上 进行 磁场 强度 的 测量 ， 而 不 能 得 到 TI 上 磁场 的 连续 分 布 H(PR)，Po € 
FF。 其 次 ， 即 使 得 到 了 T 上 磁场 的 连续 分 布 ， 对 于 其 上 磁 荷 分 布 w(P) 的 确定 涉及 到 求解 积分 方 
程 (10)， 这 通常 和 直接 求解 Neumann 外 问题 (8)，(9) 是 一 样 难 的 。 所 以 ， 有 必要 对 磁场 预测 

论 公 式 进行 离散 化 ， 以 得 到 工程 中 可 直接 使 用 的 实际 预测 公式 。 

现在 考虑 积分 方程 (10)， 将 积分 曲面 工 划 分 成 个 小 曲面 Ti (7 = 1,2,… n). SETA 
曲面 的 重心 为 Pj;，P; 处 的 外 法 线 单位 方向 向 量 记 为 n;， 并 且 PP 处 的 磁场 强度 假设 已 经 由 
传感器 测 得 ， 记 为 Hi。 假 定 小 曲面 PTj 上 的 磁 荷 面 分 布 密度 函数 为 常数 wj， 于 是 积分 方 
f? (10) E P; (i = 1,2,… ,n) 处 满足 


m 
wi = -ED f [ «0 2ER asp e i m, 
了 


ja t. TPP, 
1 cos(rpp,,ni) 
— -g Yos, |f 020299 asp mti m. (12) 
j=1 Ty : 
对 于 ?27 = 1,2,…,n, $ 
cosCrpP cos(rPPna 45 
es 2T xj TPp, (13) 
这 样 ， 积 分 方程 (10) 可 以 离散 化 为 线性 方程 组 
(A 4 Dx =b, (14) 


HP A = (aij)nx»» b = (2H;:niaxi X= M ud I 是 一 个 nxn 的 单位 阵 。 
注意 一 下 系数 矩阵 A 中 的 对 角 元 aii (i = 1,2,… ,n)， 由 (13) 的 定义 可 知 ai 是 广义 积 
但 是 根据 定理 (1) 中 的 第 三 = 条 ， 可 以 证 明 as EKA. 所 以 系数 矩阵 入 中 的 每 个 元 素 都 可 以 
明确 求 出 ， 于 是 当 det(A +I) 去 0 时 ， 可 以 得 到 x = (A 十 了 -1b。 如 果 det(A I) = 0， 则 可 
以 改变 曲面 工 的 划分 方案 ， 得 到 一 个 不 同 的 A +I， 或 者 也 可 以 利用 下 述 的 矩阵 广义 道 方法 。 
由 于 在 实际 工程 中 ， 考 虑 到 传感器 的 成 本 以 及 应 用 时 的 条 件 眼 制 ， 对 于 一 个 磁场 预 
测 问 题 ， 可 能 无 法 使 用 大 量 的 传感器 ， 而 上 述 方法 所 需 的 传感器 数量 与 小 曲面 的 数量 
是 一 样 的 ， 为 n 个 。 如 果 只 用 m (m < n) 个 传感器 测 得 小 曲面 T;、，(k = 1,2,…,m) (其 
中 {1,72;… jm} C (5,2, np) 重心 处 的 磁场 强度 了 H;,， 则 可 以 取 (14) 中 的 第 庆 行 ,构成 
一 个 新 的 线性 方程 组 


(A+i)x=b. (15) 
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此 时 的 系数 矩阵 ( 信 二 了 是 mm x n 的 ， 当 m < n 时 (A+1)-1 没 有 定义 ， 因此， 可 以 转 而 
求解 (A 十 了 的 Moore-Penrose 广 义 逆 ( 久 二 iii。 根据 矩阵 论 儿 的 结果 ， 任 意 矩 阵 的 Moore- 
Penrose 广 义 逆 矩阵 存在 且 唯一 ， 并 且 可 以 通过 奇异 值 分 解 等 算法 求 出 ， 于 是 可 以 取 x = 
(A -- 1)!b. 

在 求解 出 x = (wina 之 后 ， 根 据 (11)， 区 域 Q' 中 的 磁场 强度 的 预测 公式 变 为 


1 
H(M) = 二 1 o(P) PH dsp 
T PM 


1 n 
=r? wy {| 5 asp, Mew. (16) 
= PM 
IS rj 


A (14), (16) 8X (15), (16) 即 为 磁场 预测 的 实际 公式 。 

对 于 舰 船 磁 场 预测 的 实际 公式 ， 有 几 点 需要 说 明 的 地 方 : 

1) 注意 方程 组 (14) 的 右 端 项 b 的 表达 式 ， 可 以 发 现 ， 其 实 并 不 一 定 要 知道 在 Pi (i = 
1,2,… ,n) 处 的 磁场 值 H， 而 只 需要 知道 在 这 些 点 上 磁场 的 法 向 分 量 Hi n; 即 可 。 这 对 于 工 
程 应 用 是 有 利 的， 因为 只 需 在 相应 的 点 上 沿 着 法 向 ni; 使 用 单 分 量 传感器 测量 即 可 ， 而 不 必 使 
用 三 分 量 传感器 进行 测量 ， 可 以 有 效 降低 应 用 成 本 。 

2) 利用 传统 的 深度 换算 方法 只 能 预测 半空 间 中 的 舰 船 磁场 强度 ， 然 而 在 某 些 特殊 的 应 用 
场合 下 ， 只 预测 半空 间 上 的 磁场 强度 是 远 远 不 够 的 。 比 如 对 于 潜艇 而 言 ， 除 了 要 对 抗 下 方 空间 
中 的 水 雷 ， 同 时 还 要 对 抗 上 方 的 航空 磁 探 。 然 而 本 方法 可 以 预测 任意 形状 闭 曲面 工 之 外 任何 点 
上 的 磁场 强度 ， 所 以 可 以 看 作 是 深度 换算 方法 在 任意 方向 上 的 一 种 推广 。 

3) 对 于 传统 的 磁体 模拟 法 ， 其 中 的 主要 难点 在 于 确定 等 效 源 的 位 置 ， 这 牵涉 到 对 舰 船 等 
铁 磁 物 体 结构 的 先 验 知识 。 也 就 是 说 ， 我 们 要 事先 知道 被 研究 的 舰 船 等 铁 磁 物体 的 几何 结构 ， 
才能 较 合 理 地 确定 等 效 源 的 位 置 ， 这 是 磁体 模拟 法 能 够 准确 预测 磁场 的 一 个 必要 条 件 。 然 而 在 
一 般 的 情况 下 ， 了 和 解 铁 磁 物 体 的 几何 结构 是 很 困难 甚至 是 不 可 能 的 ， 比 如 舰 船 内 部 的 铁 磁 主 
机 、 辅 机 等 。 但 是 在 本 方法 中 ， 并 不 需要 了 解 铁 磁 物 体 的 几何 结构 ， 这 就 好 比 把 舰 船 等 铁 磁 物 
体 放 在 一 个 黑 盒 内 ， 其 在 黑 盒 外 产生 的 磁场 可 以 看 成 从 黑 盒 表面 散发 出 来 ， 于 是 舰 船 等 铁 磁 物 
体现 在 就 变 成 了 等 价 的 薄 壳 黑 盒 表 面 ， 这 样 就 把 对 于 复杂 的 铁 磁 物体 系统 的 研究 转化 为 薄 壳 黑 
盒 表 面 的 研究 ， 而 后 者 是 人 为 取 定 的 ， 具 有 简单 的 几何 形状 ， 并 且 总 是 已 知 的 。 因 此 ， 本 方法 
比 磁体 模拟 法 具有 更 大 的 适用 范围 。 


4 ”数值 模拟 


考虑 位 于 给 定 直 角 坐 标 系 中 的 两 个 旋转 椭 球 体 EBl 和 和 EE2 以 及 一 个 磁 偶 极 子 D1: GER 
体 玉 1 的 长 半 轴 al = 5m， 短 半 轴 b1 = cl = 2m， 中 心 在 坐标 原点 ， 长 轴 沿 z 轴 方向 ， 总 
WS Mg; = (20000,0,0)Am?; WERI 2 的 长 半 轴 a2 = 2.5m， 短 半 轴 好 = c2 = lm， 中 心 
在 (0,0, 1)m, KHA r 轴 方 向 , 总 磁 矩 Mas = (—2500,0,0)Am?; 磁 偶 极 子 D1 位 于 (2.5, 0.5, 0)m, 
RE Mp; = (2500,0,0)Am2 。 椭 球体 和 磁 偶 极 子 的 布置 方式 如 图 2 所 示 。 这 里 磁 和 矩 大 小 的 选择 
是 使 得 椭 球 体 用 来 模拟 舰 船 感应 磁性 ， 偶 极 子 来 模拟 舰 船 固定 磁性 。 椭 球体 和 磁 偶 极 子 在 空间 
任意 点 所 产生 的 磁场 精确 值 可 以 通过 解析 式 来 计算 叫 。 
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确定 两 条 测量 线 (单位 : 米 ) 
LIA = ((z,y,z)| 10€ z < 10, y —0, z= 4}, 
L2= {(z,y,2)| -10£z < 10, y 20, z = —4}, 
在 这 两 条 线 上 进行 磁场 的 预测 ， 并 将 其 与 精确 值 比较 。 把 曲面 工 取 作 长 方 体 
{(z,y,2)| -6xz«6, -25«y «25, — 2.5 <z x25) 


的 表面 沿 z 方 向 作 8 等 分 ， 沿 y 和 z 方 向 各 作 3 等 分 ， 如 图 3， 这 样 形成 114 个 小 曲面 ， 每 
个 小 曲面 上 的 重心 、 法 向 、 重 心 处 的 磁场 强度 都 可 以 算出 。 


图 2: 试验 设计 方案 图 3: 曲面 了 的 划分 


通过 求解 式 (14)， 再 利用 (16)， 可 以 预测 出 测量 线 地 和 Z2 上 的 磁场 值 ， 其 与 精确 值 的 比 
较 如 图 4 和 5 所 小。 实际 工程 中 常用 相对 误差 


max | 精确 值 -预测 值 | }/ max( | 精确 值 | } x 10096 
来 衡量 磁场 预测 精度 。 此 时 测量 线 Ll1 上 zx,y,z 三 个 方向 的 磁场 相对 误差 分 别 为 6.96%、12.01%、 


7.65%; L2 Lx, y, z 三 个 方向 的 磁场 相对 误差 分 别 为 5.87%、12.01%、7.08%， 可 以 满足 工程 


测量 线 Lt 测量 线 L2 
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图 4: 元 1 上 磁场 的 预测 值 与 精确 值 图 5: L2 上 磁场 的 预测 值 与 精确 值 
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接 下 来 ， 减 少 测量 点 的 数目 , 荆 的 上 下 表面 上 的 测量 点 数目 不 做 改变 ， 而 在 其 侧面 上 
取 少 量 的 测量 点 ， 取 法 如 图 6 所 示 ， 只 在 灰色 的 小 曲面 上 进行 磁场 测量 ， 总 测量 点 56 个 ， 
比 先前 的 114 个 少 了 一 半 多 。 利 用 式 (15)、(16)， 测 量 线 Ll1 和 ZL2 上 的 磁场 预测 值 与 精确 值 如 
图 7 和 8 所 示 。 此 时 测量 线 L1 上 x,y,z 三 个 方向 的 磁场 相对 误差 分 别 为 9.70%、27.77%、9.61%; 
L2 上 zw,y,z 三 个 方向 的 磁场 相对 误差 分 别 为 7.47%、27.77%、8.56%。z,z 两 个 分 量 的 误差 变 
化 不 大 ， 可 以 满足 工程 需求 ， 而 y 分 量 的 误差 变化 较 大 是 由 于 其 精确 值 接 近 0 所 引起 的 ， 因 此 
也 可 以 满足 工程 需求 。 由 此 不 难看 出 ， 虽 然 我 们 减少 了 一 半 的 测量 点 ， 但 是 预测 的 精度 还 是 相 
当 高 的 ， 所 以 利用 本 方法 ， 只 要 测量 点 的 位 置 取得 比较 合理 ， 即 使 减少 一 定 的 测量 数目 ， 其 对 
预测 精确 性 的 影响 也 是 不 大 的 。 


WEHLI 测量 线 L2 


图 7: ZL 上 磁场 的 预测 值 与 精确 值 图 8: L2 上 磁场 的 预测 值 与 精确 值 


5 ARE 


本 文 利用 位 势 理 论 ， 把 对 舰 船 等 铁 磁 物体 的 磁场 预测 问题 的 研究 转化 为 对 等 效 磁 荷 面 分 布 
的 研究 。 该 方法 具有 坚实 的 数学 理论 作为 保证 ， 适 用 范围 广 ， 并 且 十 分 有 利于 工程 实现 。 数 值 
模拟 也 表明 了 此 法 的 精确 性 。 
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A Magnetic Field Prediction Method Based on Potential Theory 
CHEN Jie!, LU Xi-wen!, ZHANG Rong?, LI Jian-wu? 
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2- No.704 Institute, CSIC, Shanghai 200031) 


Abstract: To enhance the stealth ability of naval vessels, degaussing processes should be taken. The 
magnetic field prediction is a key step in degaussing. For an arbitrary naval vessel, a surface enclosing 
it is constructed. Based on the potential theory, an integral equation relating the equivalent magnetic 
charge layer on this surface with the measurements of magnetic field sensors is developed. After 
solving the equation by discretization, the magnetic charge layer is obtained to predict the magnetic 
field produced by the vessel. This method needs no knowledge of the geometric structure of the studied 
ferromagnetic objects and can make the prediction in any direction, so it can be widely applied. The 
numerical simulations show the precision and the practical value of the method. 

Keywords: naval vessels; degaussing; magnetic field; prediction; potential theory; integral equation 


